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Abstract —  We  consider  distributed  parameter  systems  where 
the  underlying  dynamics  are  spatially  periodic  on  the  real 
line.  We  examine  the  problem  of  exponential  stability,  namely 
whether  the  semigroup  decays  exponentially  in  time.  It 
is  known  that  for  distributed  systems  the  condition  that  the 
spectrum  of  A  belongs  to  the  open  left-half  plane  is,  in 
general,  not  sufficient  for  exponential  stability.  Those  systems 
for  which  this  condition  is  sufficient  are  said  to  satisfy  the 
Spectrum  Determined  Growth  Condition  (SDGC).  In  this  work 
we  separate  A  into  a  spatially  invariant  operator  and  a  spatially 
periodic  operator.  We  find  conditions  for  the  spatially  invariant 
part  to  satisfy  the  SDGC.  We  then  show  that  the  SDGC  remains 
satisfied  under  the  addition  of  the  spatially  periodic  operator,  if 
this  operator  is  small  enough  relative  to  the  spatially  invariant 
one.  A  similar  method  is  used  to  derive  conditions  which 
guarantee  that  A  has  left-half  plane  spectrum,  and  thus  the 
system  is  exponentially  stable.  The  results  are  demonstrated 
through  simple  illustrative  examples. 

1.  Introduction 

In  all  engineering  applications,  the  temporal  stability  of 
a  system  is  of  central  importance.  In  linear  systems  theory, 
assessing  exponential  stability  is  of  particular  interest.  For 
finite-dimensional  systems  (systems  with  finite-dimensional 
state-space),  exponential  decay  of  ||e'^‘||  is  guaranteed  if  the 
spectrum  of  the  A-matrix  lies  inside  the  open  left-half  of  the 
complex  plane  (open  LHP). 

The  situation  is  much  more  complicated  in  the  case  of 
infinite-dimensional  systems.  For  example,  it  is  possible  that 
the  spectrum  of  the  yl-operator  of  such  a  system  lies  inside 
the  open  LHP,  and  yet  ||e'^‘||  actually  grows  exponentially 
[l]-[3].  In  such  cases  it  is  said  that  the  spectrum-determined 
growth  condition  is  not  satisfied  [3]. 

Yet  there  exists  quite  a  wide  range  of  infinite-dimensional 
systems  for  which  the  spectrum-determined  growth  condition 
is  satisfied.  These  include  (but  are  not  limited  to)  systems  for 
which  the  A-operator  is  sectorial  (also  known  as  an  operator 
which  generates  a  holomorphic  or  analytic  semigroup)  [4]- 
[6]  or  is  a  Reisz-spectral  operator  [7].  In  this  paper  we  focus 
on  sectorial  operators. 

Thus  to  establish  exponential  stability  of  a  system,  one 
line  of  attack  would  be  to  show  simultaneously  that  (i)  A  is 
sectorial  and  (ii)  the  spectrum  of  A  lies  in  the  open  LHP.  But 
this  still  does  not  make  the  problem  trivial.  In  fact  proving 
that  an  infinite-dimensional  operator  is  sectorial,  and  then 
finding  its  spectrum,  can  in  general  be  extremely  difficult. 
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In  this  paper  we  will  be  dealing  with  a  class  of  spatially 
periodic  systems  on  the  real  line.  These  are  systems  for 
which  the  system-operators  A,  B,  and  C  are  spatially  pe¬ 
riodic  [8]  (i.e.,  they  commute  only  with  spatial  translations 
of  size  equal  to  some  real  scalar  X  >  0  called  the  period).  We 
consider  the  yl-operator  as  the  sum  of  a  spatially  invariant 
[9]  operator  A°  and  a  spatially  periodic  operator  e  E,  where 
e  is  a  complex  scalar.  Our  aim  is  to  find  conditions  under 
which  this  system  is  exponentially  stable. 

To  show  (i)  and  (ii)  we  take  an  indirect  route.  We  first 
find  conditions  on  the  spatially  invariant  operator  A°  such 
that  (i)  and  (ii)  are  satisfied.  We  then  show  that  (i)  and  (ii) 
will  remain  satisfied  if  the  spatially  periodic  operator  E  is 
“weaker”  than  A°  in  a  sense  that  we  describe  and  if  e  is  small 
enough.  The  reason  for  this  indirect  approach  is  that  (i)  and 
(ii)  are  much  easier  to  check  for  a  spatially  invariant  operator 
than  they  are  for  a  spatially  periodic  one.  All  conditions  we 
derive  are  in  the  Fourier  domain  and  can  be  checked  point- 
wise  in  the  spatial-frequency  variable  fc  G  M. 

Our  presentation  is  organized  as  follows.  We  briefly  review 
the  frequency  representation  of  spatially  periodic  operators 
in  Section  II.  We  describe  the  problem  set  up  in  Section  III. 
Section  IV  deals  with  general  notions  of  the  spectrum  and 
sectorial  operators.  Section  V  contains  the  main  contributions 
of  the  paper  and  is  divided  into  two  parts;  the  first  part  deals 
with  condition  (i)  described  above,  and  the  second  part  with 
condition  (ii).  Conclusions  and  future  directions  are  given  in 
Section  VI.  Most  proofs  and  technical  material  are  relegated 
to  the  Appendix  at  the  end  of  the  paper. 

Notation:  We  use  fc  G  M  to  characterize  the  spatial- 
frequency  variable,  also  known  as  the  wave-number.  E{T)  is 
the  spectrum  of  the  operator  T,  and  Sp(r)  its  point  spectrum, 
and  p{T)  its  resolvent  set.  To  avoid  clutter,  we  do  not 
index  norms  on  different  function/operator  spaces.  We  use 
II  •  II  to  denote  both  function  and  (induced)  operator  norms 
on  infinite-dimensional  spaces,  and  the  Euclidean  norm  for 
finite-dimensional  vectors  and  matrices;  the  difference  will 
be  clear  from  the  context.  C“  denotes  the  open  left-half  of 
the  complex  plane,  and  j  :=  \/—l.  &  is  the  closure  of  the 
set  6  C  C.  We  may  use  the  same  notation  for  a  spatially 
invariant  operator  and  its  Fourier  symbol. 

II.  Preliminaries 

In  this  section  we  briefly  discuss  the  frequency  domain 
representation  of  spatially  periodic  operators.  For  a  detailed 
account  the  reader  is  referred  to  [8]  and  [10]. 

Let  u{k)  and  y{k)  denote  the  Fourier  transforms  of  two 
spatial  functions  u{x)  and  y{x)  respectively.  If  u  and  y  are 
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III.  Problem  Setup 


related  by  a  linear  operator  y  =  Gu,  we  have 

y{x)  =  /  G{x,x)u{x)  dx 

JR 

^4  (1) 

y{k)  =  /  G{k,  k)  u{k)  dn 

Jr 

where  G  and  G  are  kernel  functions  in  the  spatial  and  Fourier 
domain,  respectively.  It  is  shown  in  [8]  [10]  that  the  most 
general  spatially  periodic  operator  can  be  represented  in  the 
Fourier  domain  as  an  operator  with  a  kernel  function  of  the 
form 

G{k,K)  —  gi{k)  S{k  —  K  —  fll),  (2) 

where  gi{k),  for  each  k,  can  in  general  be  a  matrix.  Such 
a  kernel  function  can  be  visualized  in  Figure  1.  [8]  also 

describes  how  (2)  can  be  written  as 

ye  =  Ge  ug,  9  G 

which  for  any  given  6  has  a  (bi-infinite)  matrix  representation 

gQ(0  —  Cl)  g_i(0  —  ti)  g_2(^  — 

9l(^) 

§2(9  +  Cl)  gi(9  +  Cl)  gQ(9  +  Cl) 

In  this  setting,  Ge  is  diagonal  for  spatially  invariant  opera¬ 
tors,  and  Toeplitz  for  periodic  pure  multiplication  operators. 

Example  1:  A  =  dx  and  F{x)  =  cos(nx)  have  the 
following  representations 


Ae  =  je  +  jCln 


for  every  9  G  [0,  Cl),  respectively.  Notice  that  we  have 
dropped  the  9  subscript  in  T,  as  it  is  independent  of  this 
variable.  ■ 


Let  us  now  consider  a  system  of  the  form 

dtip{t,x)  =  Aip{t,x)  +  Bu{t,x) 

=  (g1°  -f  B°  e  F  G°)  ipitjx)  +  Bu{t,x),  (4) 
y{t,x)  =Gf{t,x), 

where  t  G  [0,  oo)  and  a;  G  ffi  with  the  following  assumptions. 
The  (possibly  unbounded)  operators  A°,  B°,  G°  are  spatially 
invariant,  and  the  bounded  operators  B,  G  are  spatially 
periodic.  F{x)  =  2Lcos{Ctx)  with  L  a  constant  matrix,  and 
e  is  a  complex  scalar.  A°,  B°,  G°  and  E  :=  B°  F  G°  are  all 
defined  on  a  dense  domain  &  C  L^(IR)  and  are  closed,  u, 
y,  and  f  are  the  spatio-temporal  input,  output,  and  state  of 
the  system,  respectively. 

Then,  as  shown  in  [8]  [10]  and  also  briefly  in  the  previous 
section,  the  system  (4)  can  be  represented  as  in  the  Fourier 
domain  as  follows. 

dtMt)  =  {A°g  +  eB°gFC°g)Mt)  +  Beug{f) 

=  -f  eSg)  ipe{t)  +  BgUgit),  (5) 

ye{t)  =Cg'll)g{f), 

9  G  [0,n),  where  Bg  and  Cg  have  the  form  of  the  operator 
in  (3)  and^ 


B°g  =  B\6„)  ,  C°g 


Sg-.=  B°gFC°g  =  '■  0  ,  (6) 

Ai(0„)  0 


with  9n  '■=  9  +  Cln,  n  G  Z,  and 

Ai{-)  :=  B°{-)LC°{-  -n),  (7) 

A_,{-)  :=  B%-)LC°{-  +n).  (8) 

We  emphasize  that  the  convention  used  in  the  representation 
of  Eg  in  (6)  is  the  same  as  that  used  in  (3).  For  example  the 
n*  row  of  Eg  is  {•••  ,  0,  Ai(6'„), 0, gI_i(6»„), 0,  •  •  •  }■ 

Remark  1:  We  note  that  taking  F{x)  to  be  a  pure  cosine 
is  not  restrictive.  The  results  obtained  here  can  be  easily 
extended  to  problems  where  F{x)  is  any  periodic  function 
with  absolutely  convergent  Fourier  series  coefficients.  ■ 
Remark  2:  The  system  (4)  can  be  considered  as  the  LFT 
(linear  fractional  transformation  [11])  of  the  spatially  peri- 

’To  avoid  clutter,  we  henceforth  drop  the  from  all  Fourier  symbols. 
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Fig.  2.  The  LFT  of  a  spatially  periodic  system  G  and  a  spatially  periodic 
multiplication  operator  eF. 


odic  system  G, 


■ 

B 

B° 

G  = 

C 

0 

0 

0 

0 

and  the  (memoryless  and  bounded)  spatially  periodic  pure 
multiplication  operator  eF{x)  =  e2Lcos{^lx),  see  Figure  2. 

■ 

IV.  Spectral  &  Stability  Analysis 

It  is  shown  in  [8]  that  for  a  general  spatially  periodic 
operator  A  we  have 

S(A)  =  U  S(^e).  (9) 

ee[o,a) 

In  the  case  where  A  is  spatially  invariant  (and  thus  Ae  = 
diag  {•  •  •  ,  Ao{0n),  •••}),  (9)  further  simplifies  to 

S(A)  =  USp(Ao(fc)).  (10) 

fceM 

Example  2:  Let  A  =  Then  Ao(A:)  = 

—  see  Figure  3  (above).  Since  Ao(-)  is  scalar, 

Sp(Ao(fc))  =  Ao{k)  for  every  k.  It  is  easy  to  see  that  Ao(  • ) 
takes  every  real  negative  value  and  thus  from  (10)  A  has 
continuous  spectrum  S(A)  =  (— oo,  0  ],  see  Figure  3  (center). 

■ 

Remark  3:  When  A  is  spatially  invariant,  a  helpful  way  to 
think  about  S(A)  in  terms  of  its  symbol  Aq  is  suggested  by 
the  previous  example.  First  plot  Sp(Ao(  • ))  '^he  ‘complex- 
plane  X  spatial-frequency’  space,  as  in  Figure  3  (above)  of 
Example  2.  Then  the  orthogonal  projection  onto  the  complex 
plane  of  this  plot  would  give  S(A).  This  can  be  considered  as 
a  geometric  interpretation  of  (10).  In  Example  2,  since  Ao(  • ) 
is  real  scalar  and  takes  only  negative  values,  this  projection 
yields  only  the  negative  real  axis.  But  in  general  if  Ao(  • )  G 
C<?X9,  this  projection  would  lead  to  q  curves  in  the  complex 
plane. 

Notice  also  that  in  this  setting,  is  the  projection 

onto  the  complex  plane  of  samples  of  Sp(Ao(-))  taken 
at  k  =  9  +  Eln,  n  G  Z,  in  the  ‘complex-plane  x  spatial- 
frequency’  space.  As  9  varies  in  [0,12),  these  projections 
trace  out  S(A)  in  the  complex  plane.  This  can  be  considered 
as  a  geometric  interpretation  of  (9).  Eigure  3  (below)  shows 
the  said  samples  in  the  ‘complex-plane  x  spatial-frequency’ 
space  for  a  scalar  A.  ■ 

We  next  introduce  a  special  subclass  of  holomorphic  (or 
analytic)  semigroups.  The  reader  is  referred  to  [4]-[6]  for 


Re 


Im 


Fig.  3.  Above:  Representation  of  the  symbol  Ao(-)  of  Example  2  in 
‘complex-plane  x  spatial-frequency’  space.  Center:  S(A)  in  the  complex 
plane.  Below:  For  spatially  invariant  A,  the  (diagonal)  elements  of  Ae  are 
samples  of  the  Fourier  symbol  Ao(  • ). 


a  detailed  discussion.  Suppose  A  is  densely  defined,  p{A) 
contains  a  sector  of  the  complex  plane  |  arg(z  — a)  |  <  f +  7, 
7  >  0,  a  G  M,  and  there  exits  some  M  >  0  such  that 

||(z/-A)~^||  <1^1  for  |arg(z-a)  I  <  ^-f7.  (11) 
\z  —  a\  I 

Then  A  generates  a  holomorphic  semigroup  and  we  write 
A  G  [6]  [4].  We  say  that  A  is  sectorial  if  A 

belongs  to  some 

Einally,  a  semigroup  is  called  exponentially  stable  if  there 
exist  positive  constants  M  and  g  such  that  [7] 

||e^‘||<Me-^‘  for  t>0. 

Theorem  1:  Assume  that  A  is  sectorial.  Then  if  S(A)  C 
C“,  A  generates  an  exponentially  stable  semigroup.  ■ 
Proof:  If  A  is  sectorial  it  defines  a  holomorphic 
semigroup,  and  thus  is  differentiable  for  f  >  0  [5].  Then 
[3]  shows  that  this  is  sufficient  for  the  spectrum-determined 
growth  condition  to  hold.  In  particular,  if  S(A)  C  C“,  A 
generates  an  exponentially  stable  semigroup.  ■ 


V.  Stability  and  the  Spectrum-Determined 
Growth  Condition 

In  the  literature  on  semigroups,  there  exist  examples  in 
which  S(A)  lies  entirely  inside  C“,  but  ||e^*||  does  not 
decay  exponetially.  See  [1]  and  more  recently  [2].  In  such 
cases  it  is  said  that  the  semigroup  does  not  satisfy  the 
spectrum-determined  growth  condition  [3].  The  determining 
factor  in  the  examples  presented  in  [1]  and  [2]  can  be 
interpreted  as  the  accumulation  of  the  eigenvalues  of  Ae 
around  ±joo  in  the  form  of  Jordan  blocks  of  ever-increasing 
size  (i.e.  as  the  eigenvalues  tend  to  ±joo  their  algebraic 
multiplicity  increases  while  their  geometric  multiplicity  stays 
equal  to  one).  But  such  cases  are  ruled  out  when  one 
deals  with  holomorphic  semigroups,  which  is  the  reason  we 
consider  these  semigroups  in  Theorem  1. 

Our  ultimate  aim  in  this  section  is  to  verify  exponential 
stability.  By  Theorem  1,  in  order  to  prove  exponential  stabil¬ 
ity  of  a  holomorphic  semigroup  with  infinitesimal  generator 
A,  it  is  sufficient  to  show  that  S(A)  C  C“.  Hence,  in  the  first 
part  of  this  section,  we  give  conditions  under  which  the  A 
operators  described  by  (4)  do  indeed  generate  holomorphic 
semigroups.  In  the  second  part,  we  find  sufficient  conditions 
which  guarantee  S(A)  C  C“. 

Once  again,  the  setup  is  that  of  (4).  In  addition  assume 
that  Ao(fc)  G  is  diagonalizable  for  every  fc  S  K. 

Conditions  for  Sectorial  A 

To  find  conditions  under  which  A  in  (4)  will  define  a 
holomorphic  semigroup  we  again  use  perturbation  theory. 
We  first  find  conditions  under  which  A°  is  sectorial.  We  then 
show  that  A  =  A°  eE  remains  sectorial  if  E  is  ‘weaker’ 
than  A°  in  a  certain  sense  we  will  describe  and  if  e  is  small 
enough. 

In  the  next  theorem  we  present  a  condition  for  a  spatially 
invariant  with  symbol  Ao{k)  to  be  sectorial. 

Theorem  2:  Let  Ao(fc)  be  diagonalizable  for  every  k  G 
K,  and  let  I7(k)  be  the  transformation  that  diagonal¬ 
izes  Ao(k),  Ao(k)  =  U(k)  A(k)  U~^(k).  Let  K(k)  := 
||[/(fc)||  ||C/“^(fc)||  denote  the  condition  number  of  U(k).  If 
suPfcgK  K(k)  <  oo,  and  for  every  k  G  R,  p(Ao(k))  contains 
a  sector  of  the  complex  plane  |  arg(z  —  a)|  <  f-l-7,  7>0 
and  a  S  M  both  independent  of  k,  then  A°  is  sectorial. 

Proof:  See  Appendix.  ■ 

This  theorem  has  a  particularly  simple  interpretation  when 
Ao(  • )  is  scalar.  In  this  case  k  (U (fc))  =  1  for  all  k  gR.  Now 
since  Ao(  • )  traces  a  curve  in  the  complex  plane,  by  Theorem 
(2)  if  this  curve  stays  outside  some  sector  |  arg(z  —  a)  |  < 
I  -f  7,  7  >  0,  of  the  complex  plane  then  A°  is  sectorial. 

The  following  theorem  is  from  [6].  It  uses  the  notion  of 
relative  boundedness  of  one  unbounded  operator  with  respect 
to  another  unbounded  operator  [4]. 

Theorems:  Suppose  G  and  E  = 

B°  E  C°  is  relatively  bounded  with  respect  to  A°  so  that 

\\Ef\\  <  am  +  b\\A°f\\,  (12) 

with  0  <  a  <  00  and  0  <  6|e|  <  1/(1  -f  M).  Then  A  = 
A°  eE  is  a  sectorial  operator.  ■ 


This  theorem  says  that  if  A°  is  sectorial,  then  so  is  A  = 
A°-\-eEif  Eis  weaker  than  A°  in  the  sense  of  (12)  and  if  |e| 
is  small  enough.  Notice  that  at  this  point,  condition  (12)  can 
not  be  reduced  to  a  condition  in  terms  of  Fourier  symbols 
(i.e.  a  condition  that  can  be  checked  pointwise  in  k)  as  in 
Theorem  2.  This  is  because  E  is  not  spatially  invariant.  But 
once  the  exact  form  of  the  operators  B°  and  C°  is  known, 
(12)  can  be  simplified  to  a  condition  on  the  symbols  of  A°, 
B°  and  C°.  Let  us  clarify  this  statement  with  the  aid  of  an 
example. 

Example  3:  Consider  the  periodic  PDF 

dtf  =  —  c  Ip  -\-  e  dx  cos{nx)  tp  -G  u 

y  =  tp. 

It  is  easy  to  see  that  A°  =  —(d^  -f  —  c,  B°  =  dx  and 
C°  =  5{x)  (the  identity  convolution  operator).  By  formal 
differentiation  we  have 

Exp  =  dx  cos{Clx)  xp  =  —  H  sin(Ha;)  1/ -f  cos(Hx)  Oa;  t/. 

Using  the  triangle  inequality  and  ||sin(Ux)||  = 

II  cos(Hx)||  =  1  we  have 

||l;i/||  <  |H|||t/||  +  115,1/11-  (13) 

Thus  we  have  effectively  ‘commuted  out’  the  (bounded) 
spatially  periodic  operator  in  E,  and  are  left  with  only 
spatially  invariant  operators  on  the  right  of  (13).  Hence,  after 
applying  a  Fourier  transfomation  to  the  right  side  of  (13),  a 
sufficient  condition  for  (12)  to  hold  is  that 

|H|  -F  |fc|  <  a  -F  6  |(fc^  —  +  c|,  k  G  ffi, 

which  can  be  shown  to  be  satisfied  for  large  enough  a  >|H| 
and  b  >  0.  m 

Remark  4:  The  above  example  makes  clear  the  notion  of 
E  being  ‘weaker’  than  A°  that  we  mentioned  at  the  beginning 
of  this  subsection.  If  in  Example  3  we  had  B°  =  and 
C°  =  df  and  u  p,  =  5,  then  E  would  contain  a  5*  order 
derivative,  whereas  the  highest  order  of  5,  in  A°  is  4.  This 
would  mean  that  (12)  can  not  be  satisfied  for  any  choice  of 
a  and  b.  m 

Conditions  for  S(A)  C  C“ 

The  final  step  in  establishing  exponential  stability  is  to 
show  that  S(A)  C  C“.  Unfortunately  it  is  in  general  difficult 
to  find  the  spectrum  of  an  infinite-dimensional  operator. 
Thus  we  proceed  as  follows.  We  consider  the  block-diagonal 
operators  Ag,  9  €  [0,  H).  This  allows  us  to  extend  Gersgorin- 
type  methods  (existing  for  finite-dimensional  matrices)  to 
find  bounds  on  the  location  of  ^(Ae),  Ag  =  A°0  +  eEg.  This 
in  turn  we  use  to  find  conditions  under  which  S(Ae)  C  C“, 
and  thus  S(A)  C  C“. 

In  locating  the  spectrum  of  a  finite-dimensional  matrix 
T  G  the  theory  of  Gersgorin  circles  [12]  provides 

us  with  a  region  of  the  complex  plane  that  contains  the 
eigenvalues  of  T.  This  region  is  composed  of  the  union  of  g 
disks,  the  centers  of  which  are  the  diagonal  elements  of  T, 
and  their  radii  depend  on  the  magnitude  of  the  off-diagonal 


elements  [12].  This  theory  has  also  been  extended  to  the  case 
of  finite-dimensional  block  matrices  (i.e.,  matrices  whose 
elements  are  themselves  matrices)  in  [13].  Next,  we  further 
extend  this  theory  to  include  bi-infinite  (block)  matrices  Ag. 

For  every  fc  G  M,  take  Sfe  to  be  the  set  of  complex 
numbers  z  that  satisfy 

—  AQ{ky)  <  |e|  ^||GI_i(A:)|| -f  ||^i(fc)||^ ,  (14) 

where  (Jinin(2F  —  ^o(^))  denotes  the  smallest  singular  value 
of  the  matrix  zl  —  Ao{k). 

Lemma  4:  For  every  9,  the  spectrum  of  Ag  =  A°q  +  eSg 
is  contained  in  the  set 

6g  =  \J 

Proof:  See  Appendix.  ■ 

Example  4:  Let  us  consider  the  periodic  PDF  [8] 

dtf  =  +  —  cip  +  e  cos{Llx)  dx  Ip  +  u 

y  =  f.  (15) 

Comparing  (15)  and  (4)  we  have 

Ao{k)  =  -{k^  -  -  c,  B°{k)  =  1,  C°{k)  =  jk, 

B{k)  =  l,  C{k)  =  l,  L=^. 

From  (7)-(8),  Ai{k)  =  ^{k  —  Ll),  A_i(fc)  =  ^{k  +  Ll),  and 
thus  ||A_i(fc)||  -I-  ||Ai(A:)||  =  ^l\k  —  n|  -f  |A:  -f  n|).  Hence 
(14)  leads  to 

a^,,{zl  -  Ao{k))  =  \zl  -  Ao{k)\  <^-^{\k-n\  +  \k  +  H|) 

H|e|  |A:|<H 

|fc||e|  \k\>n  ’ 

which  means  that  the  set  6g  is  composed  of  the  union  of 
disks  with  centers  at  Ao(6*„)  and  radii  H|-|-|0„-|-H|). 

Figure  4  (above  &  center)  show  &g  in  the  complex-plane  x 
spatial-frequency  space  and  in  C  respectively.^  ■ 

Remark  5:  The  set 

S,(M)  :=  {z  e  C  I  a„i„(z/  -  M)  <  £}  (16) 

=  {z  G  C  I  II (z/  —  M)ip\\  <  e  for  some  ||(p||  =  1} 
=  {z  G  C  I  z  G  Sp(M  -I-  Z)  for  some  ||Z||  <  £} 

is  called  the  e-pseudospectrum  of  the  matrix  M  [14].  Clearly 
Se/(M)  C  ^^(M)  if  £'  <  £,  and  S^(M)  =  Sp(M) 
for  £  =  0.  The  pseudospectrum  is  composed  of  small 
sets  around  the  eigenvalues  of  M.  For  instance  if  M  has 
simple  eigenvalues,  then  for  small  enough  values  of  £  the 
pseudospectrum  consists  of  disjoint  compact  and  convex 
neighborhoods  of  each  eigenvalue  [15].  Thus  for  every  k  G 
M,  (14)  defines  a  closed  region  of  C  that  includes  the 
eigenvalues  of  A[){k).  Moreover,  comparing  (16)  and  the 

^We  would  like  to  point  out  that  Figure  4  (above)  is  technically  incorrect; 
once  the  spatially  invariant  system  is  perturbed  by  a  spatially  periodic 
perturbation  it  is  no  longer  spatially  invariant  and  thus  can  not  be  fully 
represented  by  a  Fourier  symbol.  Hence  its  spectrum  can  no  longer  be 
demonstrated  in  the  complex-plane  x  spatial-frequency  space. 


Fig.  4.  Above:  The  regions  viewed  in  the  ‘complex-plane  x  spatial- 
frequency’  space  (the  disks  are  parallel  to  the  complex  plane).  Center: 
Yi[Ae)  is  contained  inside  the  union  of  the  regions  Below:  The  bold 

line  shows  E(A^)  and  the  dotted  region  contains  S(A),  A  =  A°  +  eE. 


definition  of  Q5/c  in  (14)  and  we  have  05^  =  Se(Ao(fc)), 

£=|£|(||A_i(A;)||  +  ||Ai(fc)||). 

We  now  employ  Lemma  4  to  determine  whether  S(A) 
resides  completely  inside  C“,  as  needed  to  assess  system 
stability. 

Take  to  be  the  closed  disk  of  radius  £  and  center  at  the 
origin,  and  Sfc  to  be  the  region  described  by  (14).  Define  the 
sum  of  sets  by  =  {z  |  z  =  Z1-I-Z2,  zi  G  ili,  zi  G  ili}. 

Also,  for  every  /c  G  M  let  Xmaxik)  represent  the  eigenvalue  of 
Ao{k)  with  the  maximum  real  part,  and  let  K{k)  be  defined 
as  in  Theorem  2. 

Theorem  5:  For  every  k,  Sfe  is  contained  inside 
-PDrik)  with 

r{k)  :=  |£|(||A_i(fc)||  +  ||Ai(A;)||)tt(fc). 

In  particular,  if  S(A°)  C  C“  and 

r(fc)  <  |Re(Aniax(fc))| +/3  (17) 

for  every  fc  G  K  and  some  /3  <  0  independent  of  k,  then 
S(A)  c  C-. 

Proof:  See  Appendix.  ■ 

Example  5:  Once  again  we  use  the  scalar  system  of 
Example  4.  K{k)  =  1  since  Ao{k)  is  scalar,  |Re(Amax(fc))  |  = 
|(fc^  —  c|,  and 

||A_i(fc)||  +  ||Ai(A;)||  =  l(|fc-H|  +  |fc  +  H|). 


Thus  condition  (17)  becomes 

\A  (^\k  -  n\  +  \k  +  n\)  <  \{k^-^Y  +  c\  +  f3. 

If  this  condition  is  satisfied  for  some  /?  <  0,  the  dotted 
region  in  Figure  4  (below)  will  remain  inside  C“  and  thus 

S(A)CC-.  ■ 

VI.  Conclusions  and  Future  Work 

In  this  paper  we  study  the  problem  of  exponential  stability 
for  a  class  of  spatially  periodic  systems.  We  do  this  by  (i) 
finding  conditions  under  which  the  A-operator  is  sectorial 
(i.e.,  generates  a  holomorphic  semigroup)  and  thus  satisfies 
the  spectrum-determined  growth  condition,  and  (ii)  deriving 
conditions  which  guarantee  that  A  has  open  LHP  spectrum. 

Future  work  in  this  direction  would  include  extending  this 
procedure  to  larger  classes  of  semigroups  which  also  satisfy 
the  spectrum-determined  growth  condition. 

VII.  Appendix 

Proof  of  Theorem  2 

It  is  shown  in  [16]  that  a  sufficient  condition  for  A°  to  be 
sectorial  is  that  p{A°)  contain  some  right  half  plane  {z  € 
C  I  Re(2;)  >  p},  and 

\\zizl  -  A°)-^\\  <  M  for  Re(z)  >  p, 
for  some  p,  >  0  and  M  >  1. 

Now  since  Ao{k)  G  has  simple  eigenvalues  for 

every  k,  there  exists  a  matrix  U{k)  such  that  A^fk)  = 
U{k)  A{k)  U~^{k)  with  A(fc)  a  diagonal  matrix  with  ele¬ 
ments  Ai(fc)  :=  Xi(Ao{k)),  i  =  1,  -  ■  ■  ,q.  Thus  we  have 

\\z{zl  -  A°)-^\\  <  sup  (||z(z/- Ao(A:))”^|l) 

<  sup(||C/(fc)||  ||[/-i(fc)|| 

\\z{zI-A{k))-^\\) 

/  \z\  \ 

^  tvmax  sup  (  ,.  r  S'  i  A  (1  til  )  ’ 
km  Vdist[z,Sp(Ao(A:))]/ 

where  Kmax  :=  supj.gj{  K{k). 

Let  us  now  choose  M'  =  (1  -f  M  >  1,  and 

consider  for  a  given  k  the  region  of  the  complex  plane  where 


dist[z,Sp  (Ao(A:))] 


>  M'. 


This  region  (which  contains  the  eigenvalues  Xi{k))  is  con¬ 
tained  inside  the  union  of  the  circles 


Nmax  I  \  n  \  \  —  ^  ’  *  —  1 )  ■  ■  ■  )  7 ; 

\z  -  Xi(k)\ 

which  are  themselves  contained  inside  the  larger  circles 

N-A.WI  < 


Notice  that  (Al)  describes  circles  whose  radii  increase  like 
|Ai(fc)|/M,  M  >  1,  as  their  centers  Xi{k)  become  distant 
from  the  origin.  Clearly  a  sufficient  condition  for  these  circles 
to  belong  to  some  open  half  plane  {2;  G  C  |  Re(z)  <  p}  for 
all  fc  G  K  and  large  enough  M  is  that  Sp(Ao(A:)),  fc  G  M, 
reside  outside  some  sector  |  arg(z  —  a)  |  <  f  +  7,  7  >  0,  of 
the  complex  plane. 

Finally,  if  the  circles  (Al)  are  contained  in  some  open 
half  plane  {z  G  C  |  Re(z)  <  p}  for  all  /c  G  K,  then  for 
Re(z)  >  p,  z  G  p{AQ{k))  and  we  have 

km  '"dist[z,  Sp(Ao(fc))]  / 
and  thus  ||z(z/  —  A°)“^||  <  M  for  Re(z)  >  p. 

Proof  of  Lemma  4 

We  use  n^vTIlAr  to  denote  the  {2N  -f  1)  x  {2N  +  1) 
truncation  of  an  operator  T  on  where  H^v  is  the  projection 
defined  by 


2tV+l  times 


where  /  is  the  q  x  q  identity  matrix.  We  form  the  finite¬ 
dimensional  matrix  nArAtellArl  „  with  pure  point  spec- 
trum.  Then  using  a  generalized  form  of  the  Gersgorin  Circle 
Theorem  [13]  for  finite-dimensional  (block)  matrices,  we 
conclude  that 

S(n7v  Ale  L  [J  iBe„  C  Q3e„ 

|n|<Ar  nm 

where  are  regions  of  C  defined  by  (14).  Since  this  holds 
for  all  W  >  0,  we  have  S(Ale)  C  ©e- 

Proof  of  Theorem  5 

If  U{k)  diagonalizes  A^fk),  A^fk)  =  U{k)  A{k)  U~^{k), 
and  K{k)  =  ||C/(A:)||  ||[/“^(fc)||  denotes  its  condition  number, 
then  from  [17]  the  pseudospectrum  of  Ao{k)  satisfies 

^p(^o(^))  +  ®e  Q  X^s{X^o{k))  C  Sp(Glo(fc))  -f 

(A2) 

for  all  £  >  0.  Thus  the  first  statement  follows  immediately 
from  (A2)  and  =  'Ss{Ao{k))  with  s  =  |e|  (|jA_i(fc)|j  -F 
||Ai(fc)||).  To  prove  the  second  statement,  let  Cf  denote  all 
complex  numbers  with  real  part  less  than  /7  G  M.  It  follows 
from  S(A°)  C  C"  that  S(Al^)  C  C"  for  every  6.  If  (17) 
holds  then 

C  Sp(Ao(0n))  +  '^r{9n)  C 


i  =  !,■■■  ,q. 


(Al) 


for  every  n  G  Z,  and  from  Lemma  4  we  have  T,{Ae)  C  &e  C 
Cf,  for  some  /3  <  /3'  <  0  and  every  0.  Thus  S(A)  C  C“. 
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